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Summary of Accomplishments: Modeling the Cz Growth of Silicon 

Modeling Approach 

During the past eighteen months the CREL group at Washington University 
has established Itself as a leader ir. mathematical modeling of CZ growth 
of silicon. The unique approach of this modeling group was to relate 
in a quantitative manner, using models based on first principles, 
crystal quality to operating conditions and geometric variables. Superior 
quality Si crystals need to be of uniform diameter, of high, axially 
and radially uniform, resistivity and defect free. The resistivity 
profiles can readily be related to Imbedded dopant and Impurity concen- 
trations, while the dislocations are a strong function of the Interface 
shape and residual stresses in the crystal. In turn the dopant dis- 
tributions, interface shape and diameter are governed by the heat balance 
at the interface. Thus, it can be . ncluded that crystal quality can 
directly be related, in a quantitative way, to the hydrodynamics, mass and 
heat transfer in the melt and to heat transfer in and from the crystal. 

These phenomena are dependent on: geometric variables which are fixed 

for a given puller (e.g., crucible and enclosure shape, heater shape, 
etc.), process variables which are governed by the selected operating 
schedule (e.g., melt depth, height of exposed crucible wall, etc.), 
manipulated variables which could oe used for control (e.g., heater power, 
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crystal and crucible rotation speed, crystal pulling rate, heater position, 
gas flow rate, etc.). The variables to be controlled are crystal diameter 
and interface shape since they directly affect dopant distribution and 
residual stresses. 


The systematic modular approach to this complex problem developed 

at CREL is shown schematically in Figure 1. The sketch of the puller is 

showing Figure 2. The finite element method is used for all calculations. 

The heat balance at the interface, which is the key element that governs 

crystal quality can be stated as: 

(Heat of Solidification) x (Rate of Solidification) » 

(Heat Loss into - (Heat Transfer from 

Crystal from the Melt to the Inter- 

Interface) face) 
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Here term II is given by the heat transfer model for the crystal while 
term III is determined by the hydrodynamic and heat transfer model for 
the melt. Clearly, terms II and III dictate the interface shape through 
angle a and achievable pulling rate Vg at the fixed diameter. Equation 
(1) allows the general problem to be approached in a modular, sequential 
fashion. The rigorous mo. el for the heat transfer in the crystal can be 
solved first and coupled with an approximate model for the melt. Melt 
models can be gradually improved and incorporated into the overall 
model as shown in Figure 1. 

The above approach has been used successfully by the CREL group in 
modeling CZ of silicon and it should be directly extrapolable to modeling 
of GaAs. 
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Major Accomplishments 

1. A rigorous model based on first principles was established to predict 

the heat transfer in the crystal and the Interface shape [ Ramachandran 
and Dudukovic, J. Crystal Growth 71 (1985): 399-408]. The model is the 

first to have considered properly both direct and reflected radiation in 
the crystal enclosure. The model established two dimensionless parameters: 
a dimensionless flux at the interface and a Biot number as the key factors 
in determining the interface shape. This is confirmed by experimental 
data. The effect of including reflected radiation on axial temperature 
profiles is shown :.n Figure 3. The effect is even more dramatic for inter- 
face shape. 

2. A coupled problem for the puller was solved based on an approximate 
conductor dominated heat transfer in the melt [Srivastava, Ramachandran 
and Dudukovic, J. Crystal Growth 73 (1985); 487-504]. Both, the Laplace- 
Young equator Tor the melt meniscus shape and Gebhart's enclosure theory 
for radiation were used. This paper provides a detailed account for the 
effect of all variables an interface shape. It is shown that reflected 
radiation is important especially in the vicinity of the interface. 

3. Radiation view factors for short crystals during necking and shoulder 
growth are reported for the first time [Srivastava, Ramachandran and 
Dudukovic, J. Crystal Growth (1S86) in print]. These are essential in 
simulating properly the critical Initial stages of growth. 

4. The detailed model developed and described under item 2 was simplified. 
Useful simple relations were developed for growth rate and interface 
shape as shown in Figure 4. The various parameters A, B, b, Cp C 2 > C^ 
can be correlated with operating conditions. The simple models can then 

be run on line. [Srivastava, Ramachandran and Dudukovic, J. Clectrochem. 
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Soc. Solid State Sci. Techn. 133 (1985): 1007-1015]. The success of the 

simple models in matching the predictions of the detail models for growth 
rate versus radius and for interface shape versus radius are shown in 
Figure 5 and 6, respectively. 

5. A novel concept of improving the controllability of the diameter and 
achieving decoupled diameter and interface shape control is introduced 
and numerically demonstrated [Sri/astava, Ramachandran and Dudukovic, 

J. Crystal Growth (1986) in print]. The concept relies on regulating 
the convecting cooling component of the crystal by using gas jets. 
Variations of gas flow rate or gas pressure can also be used. Predicted 
increases in productivity with increased jet cooling are illustrated in 
Figure 7. 

6. A preliminary steady state model for simulation of the hydrodynamics 

of the melt has been established [Dorsey, Ramachandran and Dudukovic, AlChE 
One Day Symp., St. Louis, April, 1985]. 
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Overall Schematic of Cz Modeling 
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SURFACE TEHPERATURS 
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Predicted Crystal Surface Temperatures Based on Gebhart and 
Stefan Models and Experimental Results of Williams and Reuser 
[from Ramachandran and Dudukovic, 

J. Crystal Growth (1985): 399-408] 



Simple Models for Growth Rate and Interface Shane 

Developed from the predictions of the detailed 
heat transfer model 
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I 

Reynolds A(ca^/hr) 

B(cn/hr) 

e ras 


Humber 



1 

- 

126.5 

18.6 

C. 78 

2 

100 

128.9 

18.0 

0.79 

3 

500 

132.1 

16.5 

0.83 

4 

1000 

137.9 

15.5 

0.88 

5 

1500 

142.7 

14.6 

0.91 

6 

3000 

159.4 

13.1 

0.89 

7 

6000 

182.3 

10.1 

0.98 

8 

9000 

198.6 

7.4 

mxm\ 

9 

12000 

217.6 

5.44 

1.14 
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Variation of Interface Shape with Crystal Radius 
and Jet-Cooling 




Reynolds 

No. 

C 1 

CjtcaT 1 ) 

xlO 

xlO 

• 

0.417 

0.762 

0.215 

100 

0.420 

0.755 

0.235 

500 

0.428 

0.750 

0.328 

1000 

0.422 

0.716 

0.203 

1500 

0.424 

0.700 

0.191 

3000 

0.434 

0.668 

0.200 

6000 

0.450 

0.622 

0.173 

9000 

0.472 

0.595 

0.229 

12000 

0.494 

0.578 

0.275 


Note: Caata 11,6,7,6 and 9 art plotted. 


Radius (cm) 


Detailed Model 
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Paul Alexander, Chairman 

Six presentations were made at this session on progress in high-efficiency 
solar cell research. 

The presentation by the University of California* Los Angeles (UCLA)* was 
on optimization methods for solar cell numerical models. A description was 
given on optimizing algorithms of solar cell models to simultaneously vary 
design variables such as impurity concentrations* front junction depth* back 
junction depth, and cell thickness to maximize predicted cell efficiency. The 
optimization methodology was demonstrated on Solar-Cell Analysis Program in One 
Dimension (SCAP1D), which is a computer program written by Purdue University* 
that relates a specific solar cell efficiency to a given set of input 
conditions such as cell thickness, impurity concentration* junction depth, etc. 
The UCLA optimization method* by iterating a set of input conditions from 
SCAP1D, provided a maximum solar cell efficiency related to ranges of input 
conditions that were derived from the SCAP1D inputs. UCLA presented curves of 
computer runs that demonstrated their optimization techniques. 

The University of Pennsylvania described their work in identifying, 
developing* and analyzing techniques for measuring bulk recombination rates, 
and surface recombination velocities and rates in all regions of high- 
efficiency silicon solar cells. Professor M. Wolfe described their recent 
work of improving the accuracy of their previously developed DC measurement 
system by adding blocked interference filters. The system has been further 
automated by writing software that completely samples the unknown solar cell 
regions with data of numerous recombination velocity and lifetime pairs. The 
results can be displayed in three dimensions and the best fit can be found 
numerically using the simplex minimization algorithm. They also described a 
theoretical methodology to analyze and compare existing dynamic measurement 
techniques. 

Stanford University presented their work on carrier transport and 
recombination parameters in heavily doped silicon. They presented data of 
minority carrier diffusivity in both p- and n-type heavily doped silicon 
covering a broad range of doping concentrations from 10 ^ to 10^0 atoms 
per cm^. One of the highlights of their results showed that minority 
carrier dif f usivities are higher by a factor of 2 in silicon compared to 
majority carrier dif fusivities. 

Professor R. Swanson, of Stanford University* described the construction 
of a *.2.2% efficient single-crystal silicon solar cell fabricated at Stanford. 
The cell dimensions were 3 x 5 mm and 100 pm thick with a base lifetime of 
500 ps. The cell featured light trapping between a texturized top surface and 
a reflective bottom surface* small point contact diffusions, alternating between 
n-type , .id p-type in a polka-dot pattern on the bottom surface* and a surface 
passivation on all surfaces between contact regions. A V 0L of 0.681 V* a 
J sc of 41.5 mA/cm^, and a fill factor of 0.786 was reported for the 22.2% 
efficient cell. 
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The University of Washington presented their work on investigating the 
passivation of silicon surfaces using SiN x . They used high-frequency 
capacicance voltage (CV) measurements to determine the interface state density 
at midgap (D ss ). They indicated that the SiN x /Si interface exhibits values 
of D ss larger for n-type substrates than for p-type substrates and that, 
further, after heat treatment, this difference becomes even more pronounced. 
Also, it was found that the interface state density for SiN x /p-type silicon is 
essentially independent of dopant concentration, and the D ss for SiNx/n-type 
silicon interface increases from 10 *^ to 10*3 states cm“2 eV~* as the 
phosphorus concentration increases from 10*® to 10*® cm“ J . These results 
suggest that SiN x can passivate a P + /n cell more effectively than an n + /p 
cell. The University of Washington also indicates that they are working on 
special structures that will allow measurement of both surface recombination 
velocity and interface state density. 

Westinghouse described their work on achieving higher efficiency cells by 
directing efforts toward identifying carrier loss mechanisms, designing new cell 
structures, and developing new processing techniques. Their work, using 
techniques such as deep-level transient spectroscopy (DLTS), laser-beam-induced 
current (LBIC), and transmission electron microscopy (TEM), indicated that 
dislocations in web material rather than twin planes were primarily responsible 
for limiting diffusion lengths in web. They described tneir work in hydrogen 
implantation which indicated that they can improve lifetimes and cell 
efficiencies from 19 to 120 pm, and 8 to 10. 3Z (no AR), respectively, by 
implanting hydrogen at 1500 eV and a beam current density of 2.0 mA/cm 2 . Some 
of their processing improvements included use of a double-layer AR coating 
(ZnS and MgF 2 > and an addition of an aluminum back surface reflector. Cells 
of more than 16Z efficiency have been achieved. 



